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The modified unrestricted effective-medium refractive index is defined as one that yields accurate values of a rep- 
resentative set of far-field scattering characteristics (including the scattering matrix) for an object made of randomly 
heterogeneous materials. We validate the concept of the modified unrestricted effective-medium refractive index by 
comparing numerically exact superposition T-matrix results for a spherical host randomly filled with a large num- 
ber of identical small inclusions and Lorenz-Mie results for a homogeneous spherical counterpart. A remarkable 
quantitative agreement between the superposition T-matrix and Lorenz-Mie scattering matrices over the entire 
range of scattering angles demonstrates unequivocally that the modified unrestricted effective-medium refractive 
index is a sound (albeit still phenomenological) concept provided that the size parameter of the inclusions is 
sufficiently small and their number is sufficiently large. Furthermore, it appears that in cases when the concept 
of the modified unrestricted effective-medium refractive index works, its actual value is close to that predicted 
by the Maxwell-Garnett mixing rule. © 2014 Optical Society of America 

OCIS codes: (290.0290) Scattering; (290.4020) Mie theory; (290.4210) Multiple scattering; (290.5850) Scattering, 
particles. 
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For more than a century [1,2], the notion of an effective 
refractive index of a heterogeneous material has re- 
mained highly attractive and widely used [3], Indeed, if 
it proves to be sufficiently accurate then it allows one 
to drastically simplify the calculation of scattering and 
absorption properties of small particles made of hetero- 
geneous materials and thereby bypass the potentially 
daunting task of running direct computer solvers of 
the macroscopic Maxwell equations (MMEs) [4-6] . 
Especially significant is the notion of the unrestricted 
effective-medium approximation (UEMA) defined by 
Bohren [7] “as one that yields effective dielectric func- 
tions with the same range of validity as those of media 
that are usually taken to be homogeneous (e.g., pure 
water).” Despite the potentially great importance and 
utility of the UEMA, it has remained an unproven hypoth- 
esis: neither has it been derived directly from the MMEs 
in the form of a specific recipe nor has its very feasibility 
been established numerically. The direct analytical deri- 
vation of the UEMA by spatially averaging the solution of 
the time-domain MMEs appears to be highly problematic. 
Therefore, as explained below, the main objective of this 
Letter is to take advantage of the recent generalization of 
the superposition T-matrix method (STMM) and give 
the first, to our knowledge, quantitative validation of 
the UEMA hypothesis by invoking numerically exact 
computer solutions of the frequency-domain MMEs. 

It is well known that both the time-domain MMEs and 
the derivative concept of the refractive index are an out- 
come of averaging the microscopic electromagnetic field 
over volume elements that are infinitesimally small com- 
pared to the wavelength and yet contain vast numbers of 
elementary charges [8-12], The concept of an unre- 
stricted effective refractive index (UERI) is an intuitively 
obvious extension wherein the electromagnetic field is 
averaged over volume elements that are significantly 


(but not infinitesimally) smaller than the wavelength 
and contain large (but not vast) numbers of elementary 
inhomogeneities. While the inhomogeneities are as- 
sumed to be sufficiently small, they are considered mac- 
roscopic, i.e., can be characterized by a refractive index 
that is different from that of the host material. Let us now 
assume that the resulting volume-averaged electromag- 
netic field satisfies the time-domain MMEs written for 
a homogeneous material characterized by a certain re- 
fractive index. Obviously, this artificial refractive index 
serves as the UERI of the heterogeneous material. 

It is clear that for the UEMA to work in application to a 
heterogeneous target with a fixed distribution of inclu- 
sions, each volume element used for spatial averaging 
must contain a very large number of inclusions with suf- 
ficiently random (even though fixed) positions [8-12], 
This requirement can be expected to be unnecessarily 
restrictive in many practical applications. To give the 
UEMA a better chance of success, it is appropriate to 
let the inclusions move and perform averaging over time 
t as well as over the elementary volumes. 

It should be kept in mind, however, that averaging a 
time-harmonic electromagnetic field over t yields a zero 
net result owing to the rapidly oscillating complex expo- 
nential exp(-imt), where a> is the angular frequency and 
i = (— l) 1 / 2 [13,14]: 

1 n+T 

— / d t' exp(-imt') = 0. (1) 

T J t T»2n/a> 

Therefore, one must instead perform temporal averaging 
of a representative set of optical observables defined 
such that the complex exponential exp(-imt) becomes 
canceled out upon multiplication by its complex- 
conjugate counterpart [13,14], Doing that can be further 
simplified by assuming statistical ergodicity of the 
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random heterogeneous object and replacing temporal 
averaging with ensemble averaging, i.e., by averaging 
over a statistically random and uniform distribution of 
the positions of the inclusions [13,14]. This so-called 
modified UEMA (MUEMA) is declared to work if all 
the resulting ensemble-averaged optical observables 
coincide with those calculated by solving the MMEs 
for a homogeneous object with a certain value of the re- 
fractive index. The latter quantity is called the modified 
UERI (MUERI). A conceptually similar quantity was 
called the “equivalent refractive index” in Ref. [15], 

In this Letter, we test the feasibility of the MUEMA by 
using the set of six far-field optical observables that form 
the real-valued normalized Stokes scattering matrix 
according to 
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and serve to characterize the angular distribution and 
polarization state of the scattered light in the far zone 
of a random heterogeneous object [13,14], In Eq. (2), 0 e 
[0°, 180°] is the angle between the incidence and scatter- 
ing directions, while the (1, 1) element is the conven- 
tional phase function normalized according to 

dOsinOcqC©) = 1. (3) 

The specific block-diagonal structure of the matrix (2) 
implies that the random scattering object is assumed 
to be statistically isotropic and mirror-symmetric. 

For the purposes of our analysis, we model the hetero- 
geneous object as a spherical particle randomly filled 
with N identical small spherical inclusions, as shown 
by the inset in Fig. 1(a) . The size parameter of the spheri- 
cal host X — 2nR / a is fixed at 10, where R is the particle 
radius and A is the wavelength. The radius r of the inclu- 
sions is varied such that the inclusion size parameter 
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with a fixed size parameter X = 10. The color legend in panel (b) applies to panels (a)-(f). 
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x — 2nr / a takes on values 0.3, 0.5, and 1. The refractive 
indices of the host and the inclusions are fixed at 1.33 and 
1.55, respectively. 

Densely packed inclusions separated by distances 
much shorter than the wavelength cannot be considered 
“independent” scatterers characterized by “individual” 
optical cross-sections and scattering matrices. This fac- 
tor obviously disqualifies approximate approaches such 
as the ray tracing Monte Carlo technique [16]. Therefore, 
in this Letter the elements of the scattering matrix (2) for 
the heterogeneous object are calculated using a direct 
numerically exact computer solver of the MMEs, viz., 
the recently extended STMM, which is now applicable 
to arbitrarily clustered and nested spherical domains 
[17], Unlike the numerical approach used in Ref. [15], 
the extended STMM is not based on a recursive scheme 
known to generate questionable results [ 18 , 19 ], 

Following the methodology introduced in Refs. [ 20 , 21 ], 
the statistical randomness and uniformity of the object’s 
interior is modeled in two steps. First, a random-number 
generator is used to create a fixed quasi-random and 
quasi-uniform configuration of the N inclusions while 
making sure that the volumes of the inclusions do not 
overlap and do not cross the host’s boundary. Second, 
all far-zone optical observables are averaged over the 
uniform orientation distribution of the resulting hetero- 
geneous object. 

It is obvious that if the boundary of the host particle is 
perfectly spherical then any effective-medium approxi- 
mation must satisfy the well-known Lorenz-Mie identity 
[ 22 ] 

a 2 (0)/a,(0) = l. (3) 

Therefore, a deviation of the ratio a 2 (0) / a, \ (0) from 100% 
is the most direct and unequivocal indicator of the 
numerical inaccuracy of the MUEMA. Figure 1(a) shows 
that the inclusion size parameter x — 0.3 yields 
o 2 (0) /<t \ (0) values hardly distinguishable from 100%, 
whereas the inclusion size parameter x — 1 results in 
an obvious failure of the MUEMA. The size parameter 
x — 0.5 causes deviations of the ratio a 2 (0) / a \ (0) from 
100% that are noticeable, but do not necessarily invali- 
date the MUEMA. Note that the volume density of the 
inclusions is approximately the same in all three cases: 
it is 21.6% for x — 0.3 and x — 1 and 20% for x — 0.5. 

The remaining nonzero elements of the scattering 
matrix are depicted in panels (b)-(f) of Fig. 1. The green 
curves show the results of Lorenz-Mie computations for 
a homogeneous spherical particle with X — 10 that pro- 
vide the best fit to the STMM results for a heterogeneous 
particle with x — 0.3 and N — 8000 shown by the yellow 
curves. The nearly perfect agreement provides a convinc- 
ing validation of the MUEMA hypothesis. The corre- 
sponding best-fit Lorenz-Mie refractive index value is 
TOjjy = 1.376, which is almost identical to the value 
predicted by the well-known Maxwell-Gamett effec- 
tive-medium rule [1], However, the red curves in these 
panels demonstrate that this rule definitely fails in appli- 
cation to larger inclusions with x — 1. The latter conclu- 
sion is consistent with the results reported in Ref. [20], 

In Fig. 1(g), the yellow curve shows the STMM results 
for x — 0.3 and N — 500. The best-fit Lorenz-Mie results 


are shown by the green curve and correspond to the re- 
fractive index m LM = 1.33. This best-fit value is different 
from the value TOlm = 1.3328 predicted by the Maxwell- 
Gamett mle. The red curve computed for m LM = 1.3328 
reveals a noticeably worse agreement with the STMM 
curve. This result may indicate that the number N — 
500 of inclusions may not be sufficiently large to ensure 
the requisite quasi-homogeneity of the interior of the 
heterogeneous spherical particle at spatial scales signifi- 
cantly smaller than the wavelength. 

Figure 1(h) compares the STMM results for x — 0.5 
and N = 1600 with the best-fit Lorenz-Mie results ob- 
tained for to lm = 1.375. Although this best-fit refractive 
index is close to that predicted by the Maxwell-Garnett 
rule, the deviations of the green curve from the yellow 
one are noticeable and suggest that the very concept 
of the MUEMA becomes questionable at inclusion size 
parameters comparable to and exceeding 0.5. This con- 
clusion is reinforced by Fig. l(i) which demonstrates that 
the STMM results for x — 1 and A = 216 cannot be repro- 
duced by Lorenz-Mie computations for any m LM value. 

In summary, the recently improved STMM [17] has 
finally provided a reliable direct means of testing the 
general idea of the modified UEMA. Instead of calculat- 
ing integral radiometric observables such as the optical 
cross-sections, we based our analysis on the elements of 
the Stokes scattering matrix as far-field observables 
much more sensitive to microphysical properties of 
the scattering object. Our comparisons of numerically 
exact results computed for randomly heterogeneous 
and homogeneous spherical particles have shown 
unequivocally that the MUEMA is a sound approach pro- 
vided that the size parameter of the inclusions is suffi- 
ciently small and their number is sufficiently large. 
The threshold values of x and N can be expected to de- 
pend on the refractive indices of the host and the inclu- 
sions as well as on the size parameter of the host and 
should be further analyzed and quantified. It is important, 
however, that the remarkable quantitative agreement be- 
tween the yellow and green curves in Figs. 1(b)— 1(f) over 
the entire range of scattering angles validates the general 
concept of the modified unrestricted effective refractive 
index. This Letter did not aim at deriving the MUEMA 
from the MMEs and formulating a specific prescription 
for the calculation of the MUERI based on the physical 
characteristics of the host object and the inclusions. 
Nevertheless, it appears from our limited numerical data 
that if the MUEMA is valid then the MUERI is likely to be 
close to that predicted by the Maxwell-Gamett mle. 
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